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A possible unified model of short and long gamma-ray bursts (GRBs), X-ray rich GRBs, and 
X-ray flashes is proposed. It is assumed that the jet of a GRB consists of many emitting sub-shells 
(i.e., an inhomogeneous jet model). The multiplicity of the sub-shells along a line of sight n s is an 
important parameter. If n s is large (3> 1) the event looks like a long GRB, while if n s = 1, the 
event looks like a short GRB. Finally, when n 3 = 0, the event looks like an X-ray flash or an X-ray 
rich GRB. Furthermore, our model may also explain the bimodal distributions of T90 duration of 
BATSE-GRBs. Clearly, our model predicts that short GRBs should be associated with energetic 
SNe. 



I. INTRODUCTION 

For the long gamma-ray bursts (GRBs) , the cosmo- 
logical distance, the collimated jet, the massive star 
progenitor, and the association with the supernova are 
almost established or strongly suggested [30, |3(| |6(| • 
However, for short GRBs, little is known since no af- 
terglow has been observed. Observed bimodal distri- 
bution of Tqq duration implies the origin of the short 
burst is different from that of long GRB, e.g., coales- 
cencing binary neutron stars. The origin of the X-ray 
flashes (XRFs) also remains unclear although many 
models have been proposed (see [HiJ and references 
therein). The observed event rate of short GRBs and 
XRFs are about a third of and similar to the long 
GRBs, respectively [II|22|27|. Although there may 
be a possible bias effect to these statistics, in an astro- 
physical sense, these numbers are the same or compa- 
rable. If these three phenomena arise from essentially 
different origins, the similar number of events is just 
by chance. While if these three phenomena are re- 
lated like a unified scenario for active galactic nuclei 
0- El E3 ) the similar number of events is natural and 
the ratio of the event rate tells us something about 
the geometry of the central engine. We propose a uni- 
fied model in which the central engine of short GRBs, 
long GRBs and XRFs is the same and the apparent 
differences come essentially from different viewing an- 
gles. This paper is a slight extention of ref [55| (see 
also ref |45jl. 



II. UNIFIED MODEL 

It is suggested that short GRBs are similar to the 
first 1 sec of long GRBs [3. Although short GRBs 
are harder than long GRBs (2^, this difference is 



mainly due to the difference in the low- energy spec- 
tral slope while the peak energy is similar |l2j . Fine 
temporal structures in observed light curves of short 
GRBs are similar to those of the first several seconds of 
long GRBs [3i|. Other properties, such as (V/V max ), 
the angular distribution, the energy dependence of du- 
ration and the hard-to-soft spectral evolution of short 
GRBs, are also similar to those of long GRBs [2(J. 
If short GRBs also obey the peak energy-luminosity 
relation found for the long GRBs [5^. it is suggested 
that short and long GRBs have a similar redshift dis- 
tribution |T^. 

These similarities suggest that the difference be- 
tween short and long GRBs is just the number of 
pulses, and each pulse is essentially the same [39j. As 
shown in Fig. [I] using 4Br catalogue of BATSE |35j ] . 
the fluence is roughly in proportion to the duration in 
the range of 0.01 to 1000 sec Thus, we may con- 
sider that each pulse is produced by essentially the 
same unit or the sub-jet [621, an d the GRB jet consists 
of many sub-jets. If many sub-jets point to our line 
of sight, the event looks like the long GRB while if a 
single sub-jet points to us, the event looks like a short 
GRB. Since we can observe only the angular size of 
~ 7 _1 within the GRB jet with the Lorentz factor 
7, different observers will see different number of sub- 
jets depending on the distribution of sub-jets within 
the GRB jet. Since the angular size of a causally con- 
nected region is also 7" 1 < 0.01, the opening half- 
angle of a sub-jet can be much smaller than that of 
the whole GRB jet (~ 0.1), say - 0.02. 

XRFs also appear to be related to GRBs. Softer 
and dimmer GRBs smoothly extend to the XRFs 
[iH I22I I27L |4S[ . while the peak energy- isotropic lu- 
minosity/energy relations hold for GRBs as well as 
XRFs USES. The total energy including the ra- 
dio afterglow of XRF 020903, which has a measured 
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FIG. 1: The fluence 5(50 - 300 keV) as a function of T 90 
duration for BATSE bursts from 4Br catalog (Courtesy of 
Drs. S. Michikoshi and T. Suyama). From ref [55| 



redshift, might be similar to that of GRBs |43j. Other 
properties, such as the duration, the temporal struc- 
ture and the Band spectrum of the XRFs are also 
similar to those of the GRBs, suggesting that XRFs 
are in fact soft and dim GRBs. In the sub-jet model, 
XRFs are naturally expected when our line of sight is 
off-axis to any sub-jets [II El E3, El & 



III. AN EXAMPLE OF NUMERICAL 
SIMULATION OF OUR UNIFIED MODEL 

In the following, we show a numerical simulation to 
demonstrate how an event looks so different depend- 
ing on the viewing angle in our unified model |49l l55| . 
Let us consider 7V to t = 350 sub-jets, for simplicity, 
confined in the whole GRB jet whose axis is the same 
as a i) = axis (see fig. For each sub-jet the 
emission model is the same as in E^. Let the open- 
ing half-angle of the j-th sub-jet (j = 1, • • • , iVtot) be 
A0^ b , while the opening half-angle of the whole jet 
be A#tot ■ The direction of the observer and the axis 
of the j-th sub-jet are specified by ($ bs , Vobs) and 
($ 0) , (p u} ), respectively. We assume the j-th sub-jet 
departs at time from the central engine and emits 

at radius r = r (]) and time t = t ij) = t^L + r 0) //3°'c, 
where t and r are measured in the central engine frame 
and we set t^p' = For simplicity, all sub-jets are 
assumed to have the same intrinsic properties, that 



is Af?i 



0.02 rad, 



100 



10 1 



cm, 



ag' = -1, (3f = -2.5, -jhv'f = 500 keV and the 
amplitude A {]) = const, for all j. The departure time 
of each sub-jet, i d J ep is randomly distributed between 
t = and t — idur, where tdur is the active time of 
the central engine measured in its own frame and set 
to tdur = 30 sec. The opening half-angle of the whole 
jet is set to A# to t = 0.2 rad as a typical value. We 
consider the case in which the angular distribution of 
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FIG. 2: Schematic diagram of the multiple subjet model. 
Each subjet with opening half-angle A# su b is launched 
within a cone with an opening half-angle A# to t . From 
ref IH. 



sub-jets is given by 

P(d u) ,ip U) ) dd u) dip u) oc exp[-(i? (3) /tf c )72] d§ u) dip u) 

for $ 0) < A#tot — A0 su b, where we adopt z? c = 0.1 rad 
|{58| . In this case, sub-jets are concentrated on the 
= axis (i.e., the multiplicity in the center n s ~ 
10). For our adopted parameters, sub-jets are sparsely 
distributed in the range i? c < # < A0 to t, however, the 
whole jet would be entirely filled if the sub-jets were 
uniformly distributed (i.e., the mean multiplicity n s ~ 
3). Therefore, isolated sub-jets exist near the edge of 
the whole jet with the multiplicity n s «l and there 
exists a viewing angle where no sub-jets are launched. 
Figures |3 El an d El show the angular distributions of 
sub-jets and the directions of four selected lines of 
sight, the observed time-integrated spectra, and the 
observed light curves in the X-ray and 7-ray bands, 
respectively. Note here in Figure |3 "A" represents 
the center of the whole jet and is hidden by the lines 
of sub-jets. 



A. Long GRB 

When we observe the source from the $ = axis 
(case "A"), we see spiky temporal structures (the 
upper- middle panel of Fig. 1) and E p ~ 300 keV which 
are typical for the long GRBs. We may identify case 
"A" as long GRBs. 



B. XRF and X-ray rich GRB 

When the line of sight is away from any sub-jets 
(cases "Bi" and "B 2 "), soft and dim prompt emission, 
i.e. XRFs or X-ray rich GRBs are observed with E p = 
10 ~ 20 keV and ~ 4 orders of magnitude smaller 
fluence than that of case "A" (Fig. 2). The burst 
duration is comparable to that in case "A" . These are 
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FIG. 3: The angular distribution of AT-tot = 350 sub-jets 
confined in the whole GRB jet in our simulation. The 
whole jet has the opening half-angle of A#tot = 0.2 rad. 
The sub-jets have the same intrinsic luminosity, opening 
half-angles A# su b = 0.02 rad and other properties; 7 = 
100, r = 10 14 cm, a B = -1 , Pb = -2.5 , brjv' = 500 keV. 
The axes and the angular size of sub-jets are represented 
by crosses and the dotted circles, respectively. "A" repre- 
sents the center of the whole jet and is hidden by the lines 
of sub-jets. From ref [55l| . 



quite similar to the characteristics of XRFs |19U31ll50l 
We may identify the cases "Bi" and "B 2 " as 
XRFs or X-ray rich GRBs. 

In our previous works [S3, IU 03 > we considered the 
homogeneous, instantaneous emission of the whole jet. 
Then XRFs and X-ray rich GRBs occur only when the 
whole jet is viewed off-axis (corresponding case B2). 
We now introduce sub-structure of the jet emission. 
Then, there exists the case when the observer sees 
all sub-jets off-axis but his line of sight is within the 
whole jet (case Bi). This leads two types of X-ray 
afterglows because in the afterglow phase, the viewing 
angle of the whole jet becomes important. For details 
see section 4.3. 
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FIG. 4: Time-integrated energy spectrum of the emission 
from the multiple sub-jets for the observers denoted by 
"A", "Bi", "B 2 ", and "C" in Figure 01 The source are 
located at z = 1. From ref 15511. 



may identify the case "C" as the short-hard GRBs. In 
other words, the hardness of 3 comes from as = — 1 
in our simulation so that if as ~ —0.58, the hardness 
will be 6 or so. We suggest here that not only the 
isotropic energy but also the photon index may de- 
pend on d. Another possibility is that if short GRBs 
are the first 1 sec of the activity of the central engine, 
the spectrum in the early time might be «b ~ —0.58 
for both the sub-jets in the core and the envelope. 
This is consistent with a high KS test probability for 
E p and o,b 0- These possibilities may have some- 
thing to do with the origin of as ~ —1 for the long 
GRBs. 



D. X-ray pre-cursor /post-cursor 

It is quite interesting that in Figure [5] we see the 
X-ray precursor at T Q b s ~ 60 sec in "B 2 " and the 
postcursor at T Q b s ~ 65-75 sec in "Bi". These can be 
understood by the model proposed by [3l| . 



If the line of sight is inside an isolated sub-jet 
(case "C"), its observed pulse duration is ~ 50 times 
smaller than case "A" . Contributions to the observed 
light curve from the other sub-jets are negligible, so 
that the fluence is about a hundredth of the case 
"A" . These are quite similar to the characteristics of 
short GRBs. However the hardness ratio (= ^(lOO — 
300 keV) /5(50 - 100 keV)) is about 3 which is smaller 
than the mean hardness of short GRBs (~ 6). It is 
suggested that the hardness of short GRBs is due to 
the large low-energy photon index as ~ —0.58 ^3 so 
that if the central engine launches ~ —0.58 sub- 
jets to the periphery of the core where n s is small, we 



IV. DISCUSSIONS 
A. Jet structures 

Origin of Sub-jets. — In this paper, we do not dis- 
cuss the origin of the sub-jets, but arg ue the implica- 
tions of the multiple sub-jet model |25ll3ll | This model 
is an anologue of the patchy shell model for an after- 
glow [32, Hj, • Note that in our model, a relativis- 
tically moving emitting sub-shell is called a sub-jet. 
The origin of sub-jets is not yet clear. One possibility 
is that they may arise from relativistically outflow- 
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FIG. 5: The observed X-ray and 7-ray light curves from the multiple sub-jets, corresponding the cases "A" (the upper 
left), "Bi" (the upper right), "B2"(the lower left) and "C"(the lower right) in Figure^ The sources are located at z = 1. 
From ref 15511. 



ing blobs generated by various fluid instabilities like 
Kelvin-Helmholtz, Rayleigh- Taylor instability, and so 

on0,[Hi3 

Jet structures (sub-jet configurationa) could be de- 
termined in order to reproduce all of the observed sta- 
tistical properties even including those of afterglows or 
gravitational waves |4l| ■ 

The Number of Isolated Sub-jets. — Let A# su b, $ c 
and n s be the typical opening half-angle of the sub-jet, 
the core size of the whole jet and the mean multiplic- 
ity in the core. Then the total number of the sub-jets 
(JVtot) is estimated as N tot = n s (d c / A6 suh ) 2 ~ 10 3 , 
so that the total energy of each sub-jet is ~ 10 48 erg. 
In our model, the event rate of long GRBs is in pro- 
portion to Let M be the number of sub-jets in 
the envelope of the core with a multiplicity n s = 1. 
Then the event rate of short GRBs is in proportion 
to MA#g ub , so that M ~ 10 is enough to explain the 
event rate of short GRBs. 

Angular Distribution of Sub-jets. — Of course, the 
above numerical values are typical ones and should 
have a dispersion [^J. Our core-envelope sub-jet 
model can have a similar structure to the two com- 



ponent jet model 0, 0, |H, by varying such as 
n s and M. However the distribution of sub-jets could 
also have other possibilities, e.g., a hollow-cone dis- 
tribution like a pulsar, a power law distribution, a 
Gaussian distribution 40, 58, 59] and so on. 



B. Properties of Prompt Emission 

Bimodal Distribution of T90 Duration. — It is also 
found that our model can reproduce the bimodal dis- 
tribution of T90 duration of GRBs observed by BATSE 
[Ifil ] . In our model, the duration of n s = 1 burst is de- 
termined by the angular spreading time of one sub-jet 
emission, while that of n s > 2 burst is determined 
by the time interval between the observed first pulse 
and the last one (see also Fig. 2 of ref These 
two different time scales naturally lead a division of 
the burst T90 durations into the short and long ones. 
We show an example in Figure. The dispersion of 
the lognormal-like distribution seems relatively small 
compared with the observations. This is ascribed to 
a simple modeling in this paper. We fix the jet con- 
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figuration and use the same intrinsic properties of the 
subjets. If we vary i^ur for each source and 7 <J> for 
each subjet randomly, for example, the dispersion of 
lognormal-like T go duration distribution will increase 
from the general argument that the dispersion of the 
lognormal distribution increases with the increase of 
the number of the associated random variables |2fJ . In 
more realistic modeling the observed dispersion will be 
reproduced. 

It has commonly been said that the observed bi- 
modal distribution of Tgo durations of BATSE bursts 
shows the different origins of short and long GRBs. 
However, the bimodal distribution is also available as 
a natural consequence of our unified model of short 
and long GRBs. 

Temporal Structures of Long GRBs and XRFs. — 
There are three important time scales. The first is 
the duration of the central engine measured by the 
observer Td m - = (1 + z) t^. The second is the ob- 
served pulse duration of j-th sub-jet, 5T U) , which is 
given by the angular spreading time scale of each sub- 
jet. Since emission components far from the viewing 
angle of 7 _1 is dim due to the relativistic beaming 
effect, it can be approximated as 



ST { 



a + z) 



cos 9 — cos(0 + 7 1 ) 



~ (1 + ^)^(1 + 27^) , 

where 9 = max{0, 9^ — A0^l}, and 0<f ) is the viewing 
angle of the j-th sub-jet given by cos0</ ) = n obs • n^b 
where unit vectors n obs and rx^ ' b are specified by di- 
rections (#obs,<A)bs) an d (^^V '); respectively. The 
third is the time AT, that is the difference between 
times of arrival at the observer of photons that arise 
simultaneously at the nearest and the farthest side of 
the whole jet to the line of sight, i.e., 



AT 



(1 



cos(max{0, z? obs - A0 tot }) 



(1 



COs($obs + A0 t ot) 

N T 



y 2c7 2 

. Mobs + 7A6»tot) 2 , tfobs < A0 t ot 
x < 4(7tf obs )( 7 A0 tot ), d ohs > Afltot ' 

In the following, we assume z = for simplicity. 

Let us consider the case of $ obs = (then, n s ~ 
10 2 ). Bright pulses in all the sub-jet emissions are 
observed when 0<j°(= i? w) ) < 7 -1 . Then the ob- 
served duration of the brightest pulses can be calcu- 
lated as 5T U) ~ 0.1 sec, while AT ~ 70 sec, there- 
fore, 5T (3) < T dur < AT. Brightest period lasts 
for about Tj U i - Since we assume t^l (the departure 
time of each sub-jet at the central engine) is ran- 
domly distributed, brightest spikes in the observed 
light curve are uniformly distributed in this period. 



Since the mean time-interval between bright spikes is 
about Tdur/«-s ~ 0.3 sec and is larger than the pulse 
duration ST U} , those are separated with each other. 
This period is followed by that with a duration of 
AT — Tdur ~ 40 sec in which low-flux soft events come 
from sub-jets with large 0£*>. 

On the other hand, when all sub-jets are viewed 
off-axis, i.e., n s = 0, observed time profile is greatly 
altered. In this case, 5T U) increases with the viewing 
angle and observed light curves become very smooth. 
Since we see the periphery emission, the number of 
observed (bright) pulses is small. For XRFs or X-ray 
rich GRBs with tf obs ~ A0 to t ($ob s > A0 to t), the du- 
ration of bright emission period is comparable to or 
less than ~ Td ur (ST for smallest 9 V ), while overall 
duration is given by AT. There may even be the case 
when the mean pulse interval in the observed light 
curve is smaller than the pulse duration and observed 
sub-jet emissions overlap with each other and merge 
like one or a few pulses. Therefore, variability of XRFs 
and X-ray rich GRBs is low. This is a possible expla- 
nation of observed T^-variability relation 1231. More 
quantitative analysis is under investigation |56(. 

Fine Temporal Structures of Short GRBs. — It has 
been known that a large fraction of short bursts are 
composed of a few pulses with a duration of several 
tens of msec : 33j- At present we use a very simple 
model of a subjet, i.e, instananeous emission at a cer- 
tain radius. In reality there will be a distribution of 
Lorentz factor of subjet so that the faster one will col- 
lide to the slower ones. Observed variability of short 
bursts may arise in such a more comlicated modeling 
of subjet emissions. 

Other Properties. — Some observers could see a 
cold spot with small n s in the core to have a small geo- 
metrically corrected energy even if the total energy of 
the GRBs is the same. Thus our model may be com- 
patible with the recent claim that the total kinetic en- 
ergy has smaller dispersion than the geometrically cor- 
rected 7-ray energy 0,0- The X-ray pre-/post-cursor 
is also expected if off-axis sub-jets are ejected earlier 
(for precursor) or later (for postcursor) than the main 
sub-jets [U- The viewing angle of the sub-jets may 
also cause the luminosity-lag/variability/width rela- 
tions of the GRBs including GRB 980425 [THH. 



C. Afterglow Properties 

Early Afterglows. — At an early phase, the after- 
glow variabilities may arise from the ang ular energy 
fluctuations within the GRB jet jU l34l l37| . which 
might correspond to the inhmogeneous n s . Indeed, in 
the context of the patchy shell model, the observed 
properties of variable afterglow of GRB 021004 are 
well explained if the angular size of the patches is 
about 0.02 rad js^. This size is similar to those we 
adopt in this paper. 
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FIG. 6: Tgo duration distribution in 50-300 keV of hard 
events with observed fluence ratio S(2 — 30 keV)/S(30 — 
400 keV) < 10" ' 5 . Subjet distribution is given by Gaus- 
sian form. The source redshifts are varied according to the 
cosmic star formation rate. 

Late Phase Afterglows of XRFs. — The afterglow 
could have a different behavior between the core- 
envelope sub-jet model and the uniform jet model. In 
the uniform jet model, the afterglows of XRFs should 
resemble the orphan afterglows that initially have a 
rising light curve |l4l l5l| . An orphan afterglow may 
be actually observed in XRF 030723 [uj, but the light 
curve may peak too early [Hsj . The optical afterglow 
of XRF 020903 is not observed initially (< 0.9 days) 
but may not be consistent with the orphan afterglow 
|43|. These problems could be overcome by introduc- 
ing a Gaussian tail with a high Lorentz factor around 
the uniform jet [58| because the energy redistribution 
effects may bring the rising light curve to earlier times 
[24l IHgj . Therefore, as long as the observer points 
within or slightly off-axis to the whole jet (case Bi), 
the late phase (> 1 day) properties of XRF afterglow 
may be similar to those of long GRBs. On the other 
hand, when the whole jet is viewed far from the edge 
of the jet, such that # bs ^> A0 to t (case B 2 ), XRF 
afterglows may resemble the orphan afterglow (e.g., 
see the upper-left panel of fig. 5 in ref [24(. After- 
glows for cases Bi and B 2 may correspond o b s = 5.7° 
and Oobs = 11.5° in ref [24j, respectively). Because of 
the relativistic beaming effect, case B2-like events are 
dimmer than case Bi-like events in both the prompt 
and afterglow phase (see Figure 0J, so that they may 
be rarely observed, but we believe XRF 030723 is a 
member of such a class. Recent calculation for this 
subject can be found in [l5| . 

Late Phase Afterglows of Short GRBs. — The af- 
terglow of a short GRB may be difficult to predict 
since it could resemble both the orphan and normal af- 
terglow depending on the sub-jet configuration within 
the envelope. One situation is considered in H, where 
global afterglow emission is approximated by that of 
a Gaussian structured jet superimposed on a uniform 



jet. Then the emission is dominated by the on-axis 
sub-jet at relatively early times, while it is dominated 
by the more energetic Gaussian core at later times. 
This model assumes the void around the sub-jet that 
leads to an afterglow bump. On the other hand, if the 
sub-jet structure is produced by one whole jet through 
hydrodynamical instability (§ 4.1), the afterglow of 
short GRB may be similar to those of structured jet 
model because in this case the whole jet is entirely 
filled with kinetic energy (i.e., no void around the on- 
axis sub-jet). 



D. Predictions of our Model 

A clear prediction of our unified model is that short 
GRBs should be associated with energetic SNe|63|. 
Even if the SNe are not identified with short GRBs 
due to some observational reasons we predict that the 
spatial distribution of short GRBs in host galaxies 
should be similar to that of the long GRBs. Another 
prediction is that short GRBs should have the same 
total kinetic energies as long GRBs, which might be 
confirmed by radio calorimetry. 

Interestingly our model has predicted short XRFs 
or short X-ray rich GRBs |55|. They are observed 
when isolated sub-jets are viewed slightly off-axis. 
The observed short XRF 040924 may be a kind of 
these bursts |9j. Note that the short XRFs will be 
longer than the short GRBs since the puls e duration 
grows as the viewing angle increases |19L l50| . The 
event rate of short XRFs will depend on the configu- 
ration of the sub-jets in the envelope. Further obser- 
vations are necessary to determine the envelope struc- 
ture. 



E. Comment on the NS-NS Merger Model for 
short GRBs 

Let us assume that short GRBs arises from coa- 
lescing binary neutron stars. Then the current esti- 
mate of the coalescing rate is about 10 -4 y _1 galaxy -1 
[2lT | while the event rate of long GRBs is estimated as 
~ 10 -6 y _1 galaxy -1 . If we assume that the distance 
to short GRBs is similar to long GRBs, the isotropic 
energy of GRBs should be a hundredth of long GRBs 
because the luminosity of long and short GRBs are 
similar (Figure ^) and the duration of short burst is 
typically a hundredth of that of long GRB. From the 
event rate, the opening angle of the short GRB is an 
order of magnitude smaller than the long GRBs. The 
typical energy of long GRBs is about 10 51 erg and the 
typical opening half- angle is ~ 0.1. This suggests that 
the total energy of short GRBs is about 10 47 erg, while 
~ 10 52 erg is liberated from the coalescence of binary 
neutron star. This means that the short burst is much 
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less effective and much less energetic compared to long 
GRB. 
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